Inequivalence of time and ensemble averages in ergodic systems: exponential versus 

power-law relaxation in confinement 
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Single particle tracking has become a standard tool to investigate diffusive properties, especially 
in small systems such as biological cells. Usually the resulting time series are analyzed in terms 
of time averages over individual trajectories. Here we study confined normal as well as anomalous 
diffusion modeled by fractional Brownian motion and the fractional Langevin equation, and show 
that even for such ergodic systems time-averaged quantities behave differently from their ensemble 
averaged counterparts, irrespective of how long the measurement time becomes. Knowledge of the 
exact behavior of time averages is therefore fundamental for the proper physical interpretation of 
measured time series, in particular, for extraction of the relaxation time scale from data. 

PACS numbers: 05.40.-a,02.50.-r,05.70.Ln,87.15.Vv 



I. INTRODUCTION 

Due to recent advances in single particle tracking tech- 
niques, analyses based on single trajectory averages have 
been widely employed to study diffusion in complex sys- 
tems, e.g., of large biomolecules and tracers in living 
cells Examples include the motion in the cellular 
cytoplasm of messenger RNA molecules [2|, chromoso- 
mal loci , lipid granules [1| , , and viruses [3 , telomeres 
in cell nuclei [g| , or of protein channels in the cell mem- 
brane 0- Under the assumption of ergodicity, i.e., the 
equivalence of (long) time averages (TA) with ensemble 
averages (EA) , the physical interpretation is often based 
on the time series analysis of single trajectories. For in- 
stance, particle-to-particle diffusion properties are typi- 
cally studied via TA mean squared displacements (MSD) 
of individual time series x{t), 



52(A) = 



T- A 



T-A 



^.{t + A)-x{t) dt, (1) 



where A is the lag time and T the length of the time se- 
ries. Invoking ergodicity arguments, it is tacitly assumed 
that (52(A) corresponds to the EA MSD {x^{t)) with the 
identification i ^ A, in the limit of long measurement 
times (i.e., T — >■ oo). For free normal diffusion one can 
indeed show analytically that {x'^{t)) = S^{t) = 2Kit as 
T ^ (X) At finite T, the result for (52(A) will gen- 

erally show trajectory-to-trajectory variations. However, 
a similar equivalence still holds when (52 is avera ged over 
many individual trajectories: {x^(t)) = {S'^{t)) In 
what follows we use the symbol (52 when T — >■ oo, and {S^) 
for finite T, unless specified otherwise. For anomalous 
diffusion of the form {x'^{t)) = 2Kat°' with anomalous 
diffusion constant Ka of physical dimension cm2/sec" 



and anomalous diffusion exponent a {0 < a < 2) p^ . 
the same conclusion holds if the process is described 
by fractional Brownian motionfFBM) or the fractional 
Langevin equation (FLE) [l4jl6j . 

In contrast, disagreements between TA and EA are 
not surprising for non-ergodic processes. A prominent 
example is anomalous diffusion described by continuous 
time random walks (CTRW) with diverging characteris- 
tic waiting times [17H19| : while the EA MSD scales as 



the TA MSD grows linearly 
A for free motion [ll|, . 
Ai- 



(a;2(t)) ~ t" with < k < 1, 
with the lag time, (52(A) ~ 

Under confinement one observes (52(A) ~ A^~" instead 
of the saturation plateau of the EA [ll, [l^, HH . Recently 
it was found that the TA MSD of tracers in living cells 
indeed exhibit such CTRW behavior [4, 7]. 

Here, we show that even for ergodic processes the TA 
may differ from the EA. This a priori unexpected discrep- 
ancy arises from the fact that generally dynamic variables 
are not well-defined in the TA sense, and therefore care 
is necessary when interpreting TA based on knowledge 
about the corresponding EA. We explicitly study this ef- 
fect for stochastic processes of regular Brownian, FBM, 
and FLE types, confined in an harmonic potential. Pro- 
cesses of the FBM and FLE kind are closely associated 
with the motion of tracer molecules in viscous environ- 
ments, single file diffusion, monomer motion in polymers, 
or the relative motion of aminoacids in proteins 8] . They 
have also been identified as stochastic mechanisms for the 
tracer motion in living cells and reconstituted crowding 
systems [llSli^. 

Consider first an overdamped Brownian particle in the 
harmonic potential U{x) — kx"^ /2 of stiffness k. With 
initial position x(0) the EA MSD is 

(x^(t)) = (l-e-2'=*/^)/[/3fc], (2) 

while the EA taken over the TA MSD ^ becomes 

'MKTt)) = 2 
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FIG. 1: EA and TA MSD for a Brownian particle in an har- 
monic potential. Solid lines; Eqs. ((2]), (|3]). Symbols: simu- 
lations with 7 = 1, /3 = 1, fc = 1, time step St = 0.001, and 
measurement time T = 10^. The flat curve corresponds to 
thermal initial conditions. 



Here 7 is the friction coefBcient and (3 the Boltzniann 
factor. Both quantities are identical initially, before con- 
finement effects come into play: {x'^{t)) ^ 2Kit ^ {S^{t)) 
with Ki = l/[l3j]- However, for T, A ^> cxd (with 
T - A ^ 00) the TA MSD converges to S^A) 2/[^fc], 
twice the thermal value {x'^)th = The difference 

between TA and EA is shown in Fig. [T] for both the ana- 
lytical results and simulations, with excellent agreement. 
Note the sudden dip of the TA MSD at A sa T at the 
finit e measurem ent time T, at which the limiting behav- 
ior (52 (A T)) = (x^(r)) is observed. These features 
are generic for the definition of the TA MSD ([1]) under 
confinement, compare Ref. (TBI. [20j. 

What results will be obtained for more complicated, 
non-Brownian motion? We analyze the case of anoma- 
lous diffusion governed by FBM and FLE and show that 
the entire relaxation dynamics is significantly different 
for the TA, despite the ergodic nature of these processes. 
Knowledge of the exact behavior of TA quantities is im- 
perative for the correct physical interpretation of time 
series, in particular, to extract the relaxation time. 



II. FRACTIONAL BROWNIAN MOTION 

FBM Xa{t) in an external potential U follows the 
Langevin equation 



dxa{t) 
dt 



(4) 



driven by fractional Gaussian noise £,a{t) of zero mean 
{^a{t)) = and slowly decaying, power-law autocorrela- 
tion [ty^t') [H,!!! 



(5) 



In free space, (a;„(t)) = 2Kat" [11] ■ Note the change of 
sign in Eq. ^ between antipersistent subdiffusion < 
a < 1 and persistent superdiffusion 1 < a < 2. Different 
to subdiffusive CTRW processes with diverging waiting 
time scale, FBM does not exhibit ageing effects. In fact, 
the free space propagator is the Gaussian [2^ 



Pix,t) 



1 



exp 



(6) 



whose smooth shape contrasts the pronounced cusps 
at the initial position in subdiffusive CTRW processes 
[13, [ll]. Moreover, the propagator dH) obeys a time- local 
diffusion equation with time-dependent diffusivity |25| . 
The formal solution of the FBM-Langevin equation 

m, 



Xa{t) 



(7) 



and Eq. ([5]) yield the position autocorrelation function 

+|^{e"^-|*-*^l7(a + l,fc^i) 

+e''^l*^-*il7(a + 1,^*2) 
-e'=l*^-*il7(a + l,fc|i2-ti|)} 

-1*2 -iir+^e-'=^l*^-*il 



2(a-hl)' 



kK, 



xM{a + l;a + 2; k\t2 - ti\) 



^i"+ie-^-(*i+*^)M(a-t-l;a + 2;fcii) 
2[a + 1) 

-,rP^A^'e-^('^+'^)M{a -|- 1; a + 2; kt,). (8) 

For the EA MSD we then find 
K. 



Kit)) 



-/{a + l,kt) + 2Kat°'e 



kK, 



-t"+ie-2'=*M(a -f 1; a + 2; fci), (9) 



a + l 



where 7(2, a;) = Jg dte is the incomplete 7 function 

and 



M{a;b; z) 



r(5 - a)T{a) 



e'H''-\l-t)''-''-^dt 

(10) 

is the Kummer function |221] . Fig. [2] shows simulations 
of FBM in an harmonic potential for various a val- 
ues, demonstrating excellent agreement with result ([9|). 
Asymptotic expansion of Eq. © at short times t <^k~^ 
yields free anomalous diffusion {x'^{t)) ^ 2Kat". Close 
to stationarity, we find 



Kit)) ^ {xD,^ - ^a{a - l)K^t--^e-^\ (11) 
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FIG. 2: EA MSD {xi(t)) for FBM in an harmonic potential. 
Solid lines: Eq. Symbols: simulations with parameters 
k = 0.01, T = 2048, and xo = 0. 



FIG. 3: EA and TA MSD, {xl{t)}/K^ and 5^{A)/Kc for 
FBM in an harmonic potential (a = 0.40, 1.0, 1.4, 1.60, bot- 
tom to top). Solid and dashed lines: analytical results @ 
and H13p . Symbols: simulations. Parameters as in Fig. [2] 



exponentially approaching the stationary value 



El 



Tr(« + i) 



(12) 



with the single characteristic time scale in the expo- 
nential function: as observed in Fig. beyond t > 
the stationary state is attained independent of a. This 
property enables one to study the confinement effect by 
analyzing the relaxation of {x'^{t)). We also note an in- 
teresting feature of the relaxation dynamics in the in- 
termediate timescale: somewhat counterintuitively the 
subdifFusive particle overshoots {x'^)th before a depres- 
sion back to this value, while for superdiffusion we ob- 
serve a monotonic increase (the sign of the second term 
in Eq. pT|) depends on a). Note that the a-dependence 
of the plateau value (|12p . a reminder of the fact that 
FBM is driven by an external noise and thus not sub- 
ject to the fluctuation dissipation theorem, in contrast to 
FLE motion discussed below. Phenonienologically, both 
processes are very similar. 

For the TA MSD in the limit T ^ oo we obtain the 
expression 



(52(A) = 2if„r(a + l)/fc" + 2/f„A" 



k 

kK, 



^ {e'=^r(a + 1, fcA) + e-''^r{a + 1)} 



a+1 



M{a + l;a + 2;kA), (13) 



where T{z,x) ~ dte~*t''~^ is the complementary in- 
complete 7 function. Comparison with the EA MSD, 
Eq. demonstrates a completely different functional 
behavior over all time scales, except in the short-time 
limit for which confinement is negligible. In particular, 
at A — >• oo, we find 6^{A) = 2{xl^)th for all a. 

The fundamental difference of the relaxation dynamics 
of (52(A) and {x'^{t)) is evidenced in Fig. |31 in excellent 



agreement with Eqs. (jH) and (IT3|) . In contrast to the 
exponential relaxation of Eq. ©, the TA MSD shows a 
power- law approach to the limiting value 2{x'^)th, except 
for the Brownian limit a = 1. This is manifested in the 
asymptotic form of (52(A) at A — oo. 



52(A) ^ 2{xi),,, 



Ka,T{a + 1) 



-fcA 



fc2 



2a{a - l)Ka 
fc2A2-" 



(14) 



The transient second term becomes the leading order at 
a — I. Surprisingly, in Eq. (|14p the relaxation dynamics 
is determined by the power exponent a — 2. Moreover, 
no characteristic time scale exists beyond which the MSD 
could be regarded saturated. For subdiffusion, as the al- 
gebraic decay is relatively fast A~'^ with 1 < k < 2), 
the MSD appears saturated at sufficiently long measure- 
ment time. However, the superdiffusive MSD relaxes very 
slowly as a is closer to 2 (~ A"*^ with < k < 1). Due 
to this the corresponding MSD does not show satura- 
tion even at long measurement time T. Only in the limit 
A — >■ T the TA dips back to the plateau value of the EA. 
In typical experiments, however, this feature is obscured 
by poor statistics, and thus the relaxation time would 
likely be overestimated from the TA MSD. 



III. FRACTIONAL LANGEVIN EQUATION 



The FLE describes ergodic anomalous diffusion and 
fulfills the fluctuation-dissipation theorem [l3|. In the 
potential U, the FLE motion ya{t) follows the dynamic 



4 



equation []J, l26|, [28 







at' 



+ V7/[«(a-l)/3ifa]ea ) , 



(15) 



where ^q, (t) represents fractional Gaussian noise, m is the 
particle mass and 7 the generalized friction coefHcient. In 
the FLE, the dynamic exponent of the noise is restricted 
to 1 < a < 2. This persistent noise results in subdiffusive 
motion of the FLE in the o verda mped limit. For unbiased 
motion (fc = 0), {yl{t)) = S^{t) at T ^ 00 and 



2A2 

13171 



-T[a 



(16) 



Ea.7i{z) being a generalized Mittag-LefRer function. The 
latter is defined via its Laplace image 



Ep^s (-77 t ) 



1 



(17) 



In terms of a series expansion around z = and z — )■ 00 
this function reads j27i] 



E^ 



p,5 



r{S + pn) 



E 



\ r((5 - pn) ' 



(18) 



The MSD (fTB)) accordingly turns from ballistic motion 
~ A^ at short A to subdiffusion ~ A^^" at long A 
[H [H]. In the presence of the potential, the FLE ^ 
can be solved analyticall y in the overdamped limit, the 
stationary state yielding |28| 



{ya{tl)Va{t2))th 



1 



-Ey-r 



7r(a- 1) 



1*2 -ill 



2-a 



(19) 

with E2-a{z) = i?2-Q,i- Thus, (2/^(t)) has the station- 
ary value {ya)th — ^/{f3k) for any a, contrasting the 
a-dependent result ([T^ for FBM. Moreover for T — 00 
we obtain the TA MSD 



<52(A) = 2(y^)th l-i?2-a 



k 



jr{a - 1) 



(20) 



for A ^ Tc, where the momentum relaxation time is l2f 



r(a + 3) 



2r(a - l)(2"+i - 1)7 



l/a 



(21) 



The TA MSD (gHl) behaves distinctly different from its 
EA counterpart as weh as the TA MSD (US]) for FBM: 
the TA MSD ^ grows like ~ A^"" at intermediate 
lag time, and eventually converges to 2{y'^)th for aU a 
as A — > 00. Similar to our above observations, the long- 
time behavior of Eq. (j20p exhibits a power-law relaxation, 
namely. 



52(A) « 2(y^)th (1 



7 



fcA2 



(22) 
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FIG. 4: TA MSD for FLE motion in an harmonic potential, 
for a = 1.2, 1.4, 1.6, and 1.8 (bottom to top), and T = 2". 
A representative EA MSD is included for a = 1.2. Symbols: 
Simulations with time step 5t — 0.001, stiffness k = 100, mass 
m = 1, friction coefficient 7 = 100, and P — 1, with equilib- 
rium distribution of initial position ya{0) and velocity yci(O). 
Full lines: theoretical result (|20p . Dotted lines: unbiased mo- 
tion, Eq. p6p . with momentum relaxation at a = 1.4 and 
1.8, illustrating ballistic scaling ~ at A < Tc- 



As for FBM, the dynamic exponent of the TA MSD is in- 
dependent of the confinement (k). Intriguingly the speed 
of convergence is slower as the driving noise is more 
persistent (i.e., when a 2). Therefore, opposite to the 
tendency shown in Fig. [3] for FBM, in the FLE case the 
slow particle appears not to approach 2 (y^) th- 
in Fig. m we further analyze FLE motion in an har- 
monic potential in terms of the TA MSD for various a. 
For times A ^ Tc and in the overdamped limit, the ana- 
lytical form ([20]) agrees well with the simulations results 
for all cases. While the TA MSD approaches the thermal 
value 2(?/^)th algebraically, in our simulation the slow- 
est subdiffusive case (corresponding to a = 1.8) does not 
show saturation even for long measurement times. At 
times less than Tc the TA MSD shows quadratic scaling 
^2 ^ ^2^ The dynamics within this time range is ex- 
plained well by the full solution of the free space motion 
P^. shown for a — 1.2 and 1.8. An important fea- 
ture resulting from this inertia effect are the oscillations 
of (52 that are particularly pronounced as a approaches 
2. These oscillations are intrinsic in the sense that they 
occur regardless of the confinement, due to the strong 
persistence in and inertia effects [29| . 



IV. DISTRIBUTION OF TIME AVERAGED 
MEAN SQUARED DISPLACEMENTS 



At finite sampling time T the TA MSD (52 (A, T) is 
a random variable, even for ergodic processes such as 
Brownian motion, FBM, and FLE motion. In practice, 
this means that 6^{A,T) shows pronounced trajectory- 



5 



to-trajectory variations. This stochasticity of 5"^ is mea- 
sured by the scatter probability density 0(^), in which 
the dimensionless variable is defined through [ll| 

C=F/(F). (23) 

Such scatter distributions are of Levy stable type for sub- 
diffusive CTRW processes with diverging characteristic 
time scale [ll|, [Tl, IM US IMl ■ In FigS- EE we show 
for free and confined FBM and FLE motion, with fixed 
bin size 0.1. Each graph shows the distributions at three 
different lag times A and a for given measurement time 
T. In each graph two sets of curves were shifted upwards 
for comparison, the shift value is indicated in the graphs. 



FLE were chosen as the corresponding equilibrium distri- 
butions. For the case of long measurement time (Right), 
the distributions again appear insensitive to a and A for 
the given bin size. 

C. FBM & FLE motion in free space 

To appreciate the effect of confinement on the distri- 
bution, we simulated free FBM and FLE motion for the 
same parameters as in Figs. [3] and 21 As shown in Fig. [7| 
in both cases, the scatter distributions manifest features 
consistent with the confined cases. Only the width of 
the distribution becomes narrower by the presence of the 
confining potential, in particular, at longer lag times. 



A. FBM in an harmonic potential 



D. Influence of bin size 



On the Left of Fig. [5]we present the distributions of 5'^ 
for the data shown in Fig. [31 As expected from our pre- 
vious study [13], the distributions are centered around 
the ergodic value £, — I and become wider as the lag 
time A increases. The wider distribution at longer lag 
time means that the single TA MSD trajectories tend to 
be more erratic as A approaches T. A new finding is 
that at a fixed lag time A > 1 the scatter distribution 
becomes broader as the motion is faster (i.e., growing 
a). This behavior is mainly due to the inherent prop- 
erty of FBM itself, as the same tendency is also found 
without potential (see left panel in Fig. [7]). This depen- 
dence on a is attributed to the fact that FBM is a Gaus- 
sian stationary process in which the spatial displacement 
X for time difference A is governed by the distribution 
~ exp(— a;^/[4i4r//A"]). On the Right we show the cor- 
responding distributions when the measurement time is 
increased to T = 2^^. The distributions now appear in- 
sensitive to a and A. The fact that they are less sharp 
than the analytically predicted Gaussian is due to the 
finite size effect of the binning, see below. 



B. FLE motion in an harmonic potential 

The left and right panels of Fig. [51 depict the distribu- 
tions corresponding to the parameters used in Fig. 21 with 
T = 2^^ and 2^^, respectively. Note that the variation of 
a was restricted to the range [1.0, 2.0], as FLE motion is 
only well defined for subdiffusion. The anomalous diffu- 
sion exponent in this case is given by k = 2 — a in terms of 
the scaling exponent a of the fractional Gaussian noise. 
Generally the distributions are bell-shaped. Notably, at 
fixed lag time the distribution of FLE motion tends to 
be wider as the overdamped motion becomes slower (i.e., 
for increasing a), as opposed to the case of FBM. We 
also observe that the distributions of FLE motion appear 
generally wider than those of FBM (e.g.. Fig. [51 Left), for 
the fact that the initial values of position and velocity for 



Typically experimental probing windows are limited, 
and meaningful quantitative analysis requires more or 
less coarse binning. In the context of the scatter plots 
shown here this practically means that for a bin size of 
0.1 the peak cannot exceed the value 10, for reasons of 
normalization, <!>{£,) — 1- For simulations data we 
can arbitrarily increase the accuracy and thus reduce bin 
sizes while still maintaining good statistics. Such a result 
is shown in Fig. [51 for bin size 0.01. For this resolution 
we may compare the shape of the distribution de- 
termined from simulations with the theoretical approxi- 
mation valid for short lag times A, 



^-^exp(-«-"y-^'|, (24) 



Attt* 



At* 



as derived in Ref. 3(|. Here, the scale r* is only in- 
troduced to account for correct dimensions and can be 
taken to one [time unit], compare Ref. [30j . The graph 
shows nice agreement with the measured data from the 
simulations. Interestingly the agreement is somewhat 
better in the confined case. At larger A the distribu- 
tion becomes wider than predicted by Eq. ((24)) . due to 
the strong correlation in successive square displacements, 
[x{t + A) — x{t)Y' (see Ref. [13]), and the curves split up 
for the different a. 



V. DISCUSSION 

Studying the representative example of ergodic FBM 
and FLE motion in an external harmonic potential we 
demonstrated that the TA MSD behaves significantly dif- 
ferent from the EA MSD. Thus, naive interpretation of 
single trajectory time averages based on the knowledge 
of the ensemble behavior may lead to false conclusions on 
the physics underlying the observed motion. This so far 
overlooked discrepancy is particularly relevant for the re- 
laxation behavior: while for the EA MSD the relaxation 
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5 I 

FIG. 5: Scatter distributions for fractional Brownian motion in an harmonic potential. Left: T = 2^^ = 2,048. Right: 
T = 2^^ = 13, 1072. In each graph we compare the results for three different lag times. The two upper sets of curves were 
shifted by the indicated amount, for clarity. 




0.5 1 1.5 2 0.5 1 1.5 2 



FIG. 6: Scatter distributions for fractional Langevin equation motion in an harmonic potential. Left: T — 2^^. Right: T — 2^^. 
In each graph we compare the results for three different lag times. The two upper sets of curves were shifted by the indicated 
amount, for clarity. 




5 I 

FIG. 7: Scatter distributions for free fractional Brownian and Langevin equation motion. Left: T = 2". Right: T = 2^^. In 
each graph we compare the results for three different lag times. The two upper sets of curves were shifted by the indicated 
amount, for clarity. 
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5 I 

FIG. 8: Scatter distribution for confined (Left) and free (Right) FBM with T — 2^^ , and A = 1. Data are the same as shown 
in Figs. [5]and[7]to the Left, but analyzed with bin size 0.01. As shown by the grey line the measured scatter nicely agrees with 
the predicted Gaussian distribution ((24)) . 



time can be read off directly, the corresponding TA MSD 
appears to suggest a scale-free behavior. Hence it is im- 
perative to compare to analytical or simulations results 
for the TA of the system. We note that while here we 
focused on an harmonic external potential, the findings 
reported here also pertain to other forms of confinement. 

What is the reason for this disagreement between EA 
and TA? We find that for stochastic processes converging 
to a stationary state, as {[x{t+lS.)—x{t)]'^)t\i only depends 
on A, the definition of the TA MSD ^ in the limit of 
long-time measurement leads to the general relation 



be a straightforward idea to consider 



lim (52(A,T) 



2(x2)th[l-C,(A)], 



(25) 



which is independent of diffusion models and details of 
confinement. Here 



{x{t)x{t + A))ti, 



th 



(26) 



is the normalized position autocorrelation function. 
Therefore, the time-averaged variable S'^{A) in fact is an 
indicator of the correlation of the spatial displacement, 
not the TA MSD of a trajectory. For ergodic systems 
satisfying the Khinchin theorem {Cx{A — cxd) = 0), the 

TA MSD always saturates to 6'^ — 2{x'^)th, where the 
relaxation dynamics reflecting the spatial correlation can 
be very slow although the system is already fully relaxed, 
as shown in this study. Accordingly, in performing single 
trajectory analysis, one should be aware of the poten- 
tial pitfalls in u sing ^^(A). In contrast to non-ergodic 
systems pTl [l5l . |20| [. the anomalous diffusion exponent 
a and the anomalous diffusion constant Ka can be esti- 
mated from the log-log plot of the TA MSD at short lag 
times. Meanwhile, physical quantities associated with 
confinement such as the effective confinement size and 
the relaxation time could be incorrectly deduced from 
the long- A behavior of S^{A). 

What alternative definitions of the TA MSD could be 
used to mend this problem? Instead of Eq. ^ it would 



a;2(i) 



(27) 



For ergodic processes with a stationary state, x^(t) at 
t — > oo equals 

f x'e-^^^'-^Ux/ /'e-'^^(-)dx. (28) 



However, the time dependence of x'^{t) is different from 
that of (a;2(i)). Even for free diffusion exhibiting 
(a;2(t)) = 2i^Qp, the ensemble mean of x'^(t) is 
and thus this definition does not even work for the Brow- 
nian case. If a dynamic variable like the MSD as function 
of time is concerned, it appears that no systematic way 
exists for defining a TA expectation compatible with the 
analogous EA. 

For finite measurement time T the TA MSD shows 
trajectory-to-trajectory variations, even for Brownian 
motion. Consistent with previous findings (soj for er- 
godic processes, the distributions are centered around 
the ergodic value ^ = 1. Importantly, in all cases the dis- 
tribution is almost independent of confinement, except 
for some narrowing at long lag times. For both FBM 
and FLE at long T the distributions converge, while for 
shorter T a dependence on a prevails. 

Concluding, the study of single trajectory averages is a 
non-trivial extension of the theory of stochastic processes 
knowledge of which is necessary to establish quantitative 
models for diffusion-limited processes in small complex 
systems. The current work contributes to the develop- 
ment of such a theory, and to a toolbox of diagnosis 
methods of the exact stochastic mechanism underl yin g 
experimental single particle trajectories [1, 0, [IS I32l434| . 
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